A backward Monte-Carlo method for solving parabolic partial differential equations
نویسنده
چکیده
A new Monte-Carlo method for solving linear parabolic partial differential equations is presented. Since, in this new scheme, the particles are followed backward in time, it provides great flexibility in choosing critical points in phase-space at which to concentrate the launching of particles and thereby minimizing the statistical noise of the sought solution. The trajectory of a particle, Xi(t), is given by the numerical solution to the stochastic differential equation naturally associated with the parabolic equation. The weight of a particle is given by the initial condition of the parabolic equation at the point Xi(0). Another unique advantage of this new Monte-Carlo method is that it produces a smooth solution, i.e. without δ-functions, by summing up the weights according to the Feynman-Kac formula. PACS: 02.70.Lq; 02.60.Lj
منابع مشابه
Solving parabolic stochastic partial differential equations via averaging over characteristics
The method of characteristics (the averaging over the characteristic formula) and the weak-sense numerical integration of ordinary stochastic differential equations together with the Monte Carlo technique are used to propose numerical methods for linear stochastic partial differential equations (SPDEs). Their orders of convergence in the mean-square sense and in the sense of almost sure converg...
متن کاملFeynman-Kac representation of fully nonlinear PDEs and applications
The classical Feynman-Kac formula states the connection between linear parabolic partial differential equations (PDEs), like the heat equation, and expectation of stochastic processes driven by Brownian motion. It gives then a method for solving linear PDEs by Monte Carlo simulations of random processes. The extension to (fully)nonlinear PDEs led in the recent years to important developments in...
متن کاملA Multi-fidelity Stochastic Collocation Method for Parabolic Partial Differential Equations with Random Input Data
Over the last few years there have been dramatic advances in the area of uncertainty quantification. In particular, we have seen a surge of interest in developing efficient, scalable, stable, and convergent computational methods for solving differential equations with random inputs. Stochastic collocation (SC) methods, which inherit both the ease of implementation of sampling methods like Monte...
متن کاملOn Monte Carlo algorithms applied to Dirichlet problems for parabolic operators in the setting of time-dependent domains
Dirichlet problems for second order parabolic operators in space-time domains Ω ⊂ Rn+1 are of paramount importance in analysis, partial differential equations and applied mathematics. These problems can be approached in many different ways using techniques from partial differential equations, potential theory, stochastic differential equations, stopped diffusions and Monte Carlo methods. The pe...
متن کاملMultilevel Monte Carlo method with applications to stochastic partial differential equations
In this work, the approximation of Hilbert-space-valued random variables is combined with the approximation of the expectation by a multilevel Monte Carlo method. The number of samples on the different levels of the multilevel approximation are chosen such that the errors are balanced. The overall work then decreases in the optimal case to O(h−2) if h is the error of the approximation. The mult...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2000